The phase boundary of superconducting niobium thin films with antidot arrays 
fabricated with microsphere photohthography 
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The experimental investigation of the Ic{B)-Tc{B) phase boundary of superconducting niobium 
films with large area quasihexagonal hole arrays is reported. The hole arrays were patterned with 
microsphere photolithography. We investigate the perforated niobium films by means of electrical 
directed current transport measurements close to the transition temperature Tc in perpendicularly 
applied magnetic fields. We find pronounced modulations of the critcal current with applied mag- 
netic field, which we interpret as a consequence of commensurable states between the Abrikosov 
vortex lattice and the quasihexagonal pinning array. Furthermore, we observe Little-Parks oscilla- 
tions in the critical temperature vs magnetic field. 
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Nowadays superconducting thin fihTis are used for a 
huge variety of superconducting microelectronic devices 
such as Josephson junctions, superconducting quantum 
interference devices (SQUIDs) and coplanar waveguide 
resonators. Typically, these thin films are made of type- 
II superconductors and are penetrated by quantized mag- 
netic flux when operated in magnetic fields or when bi- 
ased with sufficiently high currents. The investigation of 
these magnetic flux lines (Abrikosov vortices) and their 
individual and collective interactions with natural and 
artificial defects in the superconductor is of high interest 
and subject to many experimental and theoretical studies 
for several decades now. 

One reason for this sustained scientific attention is that 
unpinned Abrikosov vortices respond with a dissipative 
motion to any current flowing in their vicinity. In many 
cases this motion is directly related to a reduction of 
the performance (increased noise, lowered quality factor, 
shortened coherence time) of the microelectronic device. 
Defects, however, act as local energy minima and pinning 
sites and are able to reduce or even completely suppress 
vortex motion and the related dissipation [l|-|5|. For in- 
stance, it has been demonstrated that the flux noise in 
SQUIDs and the dissipation in coplanar microwave res- 
onators can be reduced by strategically positioned micro- 
holes (antidots) d 0] . 

A second and more fundamental point is associated 
with the fact, that an ensemble of Abrikosov vortices 
interacting with an ensemble of defects in the supercon- 
ductor constitutes a highly designable and controllable 
model system for repulsively interacting particles in a 
twodimensional potential landscape. In such systems it 
is possible to investigate static effects such as the for- 
mation of quasicrystals [8l-[T2| or the controlled introduc- 
tion of potential landscape disorder [l3l - [l5j as well as 



dynamic effects such as mode locking phenomena jl6l - ll8l | 
and ratchet dynamics [l^ [l^] . 

Of particular interest in both research branches is the 
case, when the typical length scales of the defect topol- 
ogy, i.e. size and mutual distance, are comparable to the 
intrinsic length scales of the superconductor, that are the 
coherence length ^ and the magnetic penetration depth 
A, which are both temperature dependent and diverge 
at the critical temperature Tc- Well below the critical 
temperature, A and ^ can usually be found in the micro- 
to nanometer range. To pattern large areas of supercon- 
ducting fllms with submicron-scaled high density arrays 
of defects constitutes a non-trivial challenge to standard 
optical (limited by resolution) or electron beam (limited 
by time) lithography. 
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FIG. 1: (Color online) (a) Layout of a 800 x 800 /im^ large 
cross shaped bridge structure with a square center area of 
200 X 200 jj.m^ for the four-probe current voltage characteri- 
zation of superconducting thin films with pinning landscapes; 
(b) zoom into (a) , niobium film with a microsphere patterned 
quasihexagonal array of antidots. 

It has been demonstrated with different approaches. 



FIG. 2: (Color online) (a) Critical current 7c vs temperature T and magnetic flux density B of a superconducting thin film 
with a quasihexagonal lattice of antidots for a voltage criterion of Vc = 5 /xV; (b) Resistance R of the niobium film vs T and B 
measured with a current / = 50 /lA. Flux density axes are normalized to the first matching flux density Bi ~ 3.4 mT. 



that taking advantage of self-assembling structures can 
provide a way out of the difficulties to cover large ar- 
eas with tiny structures on reasonable timescales. The 
techniques involved vary from using anodized aluminum 
as substrate material 21| over depositing the supercon- 
ducting film on a layer of microspheres |22| to generating 
structures by inverse diblock copolymer micelle forma- 
tion [2^. These fabrication techniques are limited to 
certain substrate materials or they induce changes in the 
substrate properties and/or the properties of the super- 
conductor. 

Here we adopt another method to fabricate large area 
quasihexagonal arrays of submicron sized antidots, which 
is independent of the substrate material and docs not 
influence the superconducting material more than any 
standard lithography process [13]. In a previous study 
we have demonstrated, that with this fabrication tech- 
nique it is possible to reduce the vortex associated losses 
in superconducting microwave resonators by more than 
one order of magnitude j25| . In the present manuscript 
we analyze the properties of our microsphere patterned 
Nb thin films by means of transport measurements close 
to the transition temperature with a particular focus on 
signatures of commensurabilities between the antidot and 
vortex lattices. We also investigate a transition between 
the wire network (width of the superconducting material 
between the antidots W < £,{T)) and the thin film regime 
with vortices {W > ^{T)) in our samples. 

We fabricated our samples by first depositing a t = 
150 nm thick niobium film on a r-cut sapphire wafer 
by dc magnetron sputtering. Afterwards we cut the 
wafer into individual chips and carried out the lithog- 
raphy steps. For the fabrication of perforated samples 



the chips were covered with photoresist and on top of 
that with a monolayer of water suspended polystyrene 
colloids in a Langmuir-Blodgett deposition process. The 
microspheres have a diameter of Dg = 770(±25) nm 
and act as a self-assembled array of UV-light focusing 
microlenses, leading to a quasihexagonal hole array af- 
ter the exposure, their removal and the resist develop- 
ment. For a perfect, hexagonal close-packed array with 
Dg = 770 nm one would get a corresponding hole density 
of Uh ~ 1.95 /xm~^. In reality, however, one could ex- 
pect deviations from the ideal packaging due to disorder 
during the self-assembling. 

After transferring the hole array into the Nb film via re- 
active ion etching (SFg) we patterned cross-shaped bridge 
structures for electric transport characterizations into the 
films. For this we used standard optical shadow-mask 
lithography and another SFg reactive ion etching step. 
Figure [T] (a) shows one of the bridge structures with a 
square center area of 200 x 200 /im^. The antidots have 
an approximate diameter Da = 370 nm, which in prin- 
ciple can be easily varied by adjusting the lithography 
exposure time. In Fig. [T] (b) a zoom-in to the niobium 
film with antidots is depicted, which shows a domain-like 
pattern of holes with some blemishes. 

Due to the cross shape of our bridges, the transport 
current in the center area is not homogeneous but some- 
what spreadened. For ratchet devices, it has been dis- 
cussed that such a sa mple g eometry may strongly affect 
experimental results j26l - l28l |. Also in our case, this ge- 
ometry might have an effect on the absolute values of 
the measured quantities, such as the critical current Ic, 
what e.g. impedes a precise determination of critical cur- 
rent densities. However, the results presented in this 



manuscript do neither sensitively depend on the abso- 
lute values of the measured quantities and on the local 
direction of the transport current. 

We have also patterned and characterized bridges with 
100 X 100 ^m^ and 50 x 50 /xm^ large center squares, but 
the experimental results showed no dependence on the 
bridge size. Note, that the niobium chips were taken 
from the same wafer as the chips for our previous study 
on resonators (25l |. 

To characterize our samples we mount them into a 
low-temperature setup, that provides a temperature sta- 
bility AT < 1 mK, and contact them electrically with 
wire bonds. We apply a magnetic field perpendicular to 
the film plane using a superconducting coil, and monitor 
current voltage characteristics (IVCs) for many values of 
magnetic field and temperature. After collecting all IVCs 
we extract the desired information as the critical current 
Ic, the critical temperature Tc or the resistance R vs 
magnetic fiux density B and temperature T . We choose 
the threshold voltage Vc defining 1^ and the measurement 
current / for the resistance R during the evaluation. To 
reduce the voltage noise, we take several IVCs at each 
value for B and T and post-process the raw data (aver- 
aging and smoothing), such that we are more sensitive to 
modulations of the Ic^Tc phase boundary with respect to 
the applied magnetic flux. 

As an overview of the whole phase boundary Fig. [2] de- 
picts (a) the critical current Ic and (b) the resistance R vs 
magnetic flux density and temperature of one of our sam- 
ples close to Tc- Obviously, the phase boundary, i.e. the 
critical current IdB), the critical temperature Tc{B) and 
the resistance R{B), boundary is non- monotonous. All 
of the quantities strongly modulate with the applied fiux 
density. There are several "canyons" and "ridges" indi- 
cating commensurate states between the fiux line lattice 
and the hole array. 

For a more detailed view of the position of the max- 
ima and minima in the 3D phase boundary, it is con- 
venient to extract single data slices. Figure [3] shows 
several individual curves for the critical current Ic{B) 
corresponding to vertical cuts for constant temperatures 
through Fig. [2] (a). The fiux density axis in Fig. [2] and 
Fig. [3] is normalized to the first pronounced maximum 
in the critical current Bi « 3.4 mT. Under the assump- 
tion, that Bi corresponds to equal densities of vortices 
and holes in the sample, we find an antidot density of 
Ua = B/4>o w 1.65 /im~^, which is somewhat smaller 
than the previously calculated n/,, « 1.95 /im^^ for a lat- 
tice without any defects. This difference is most likely 
due to defects and dislocations in the pinning array, cf. 
also Fig. [11(b). 

We of course can not be sure, that Bi indeed corre- 
sponds to equal vortex and antidot densities, as previous 
theoretical and experimental studies on randomly diluted 
and disordered triangular antidot arrays have shown a 
certain variability of the matching fields [l3l - [T5l |. The 
same studies however suggest, that the amplitudes of 
the maxima get smaller and the peaks smear out, when 
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FIG. 3: (Color online) Critical current Ic vs magnetic flux 
density _B of a superconducting thin film with a quasihexago- 
nal lattice of antidots; diflferent curves correspond to difi'erent 
temperatures T. The flux density axis is normalized to the 
flrst matching flux density Bi ~ 3.4 mT. 



the dilution/disorder in the pinning array is increased. 
The observation of quite strong and sharp maxima in 
the phase boundary of our samples at B ^ Bi and also 
a.t B = B2 = 2Bi thus suggest, that we have a rather 
ordered lattice in accordance with the impression from 
optical images (cf. Fig[T|). 

We believe, that the data of Fig. [3] also show finger- 
prints of the disorder in the lattice for higher field val- 
ues B > B2. For a perfect triangular lattice, one would 
expect maxima of the critical current predominantly at 
integer multiples of Bi. However, for higher values we 
find two shoulders around B sa 3.5i?i and B w 5Bi 
and none at _B = 3i?i or _B = 4i?i. This observation 
might refiect the necessity of additional vortices at inter- 
stiti al p ositions to stabilize the vortex lattice (cf. also 
[13I Il4l|). leading to robust configurations at i? w 3.5i?i 
and B ^5Bi. 

The effect of shifted and missing peaks is related to the 
widely discussed phenomenon of pinning interstitial vor- 
tices by "caging" them between regularly pinned ones. 
The presence of interstitial vortices can induce an in- 
creased maximum number of vortices, which can occupy 
the pinning sites. This in turn can cause an increase of 
the collective pinning strength by adding vortices. These 
effects have been observed in theoretical and experimen- 
tal studies on periodic as well as quasiperiodic pinning 
arrays [l^ [29l - [33| . In all of these cases, the caged vor- 
tices have a spatial distribution with the same symmetry 
as the underlying pinning lattice. In the present study, 
however, we assume the interstitial vortices to be filling 
up the imperfect antidot lattice at the positions of miss- 
ing antidots, what is very similar to the situation in ran- 
domly diluted pinning arrays 0, [l3l- Thus, the fiux val- 
ues, at which stable vortex configurations and matching 
features in the phase boundary appear, are more related 
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to the completeness and quality of the pinning array than 
to the intrinsic symmetry of the domain-like array parts. 
Finally, although not fully identical, all these mechanisms 
are closely related to each other. 

Besides analyzing the phase boundary by taking hor- 
izontal slices, we can also take vertical slices for chosen 
currents or voltages of the 3D boundary and end up with 
the critical temperature vs magnetic field plots (or sec- 
ond critical magnetic field vs temperature, respectively) 
for different Tc {Bc2) criteria. Figure |4] shows a plot of 
Tc{B)/Tc{0) = Tc/Tco of a perforated sample (symbols) 
for a resistance criterion Rc/Rn ~ 0.5 and measured with 
an applied current / = 50 /iA. i?„ denotes the normal 
state resistance at T = 10 K. Oscillations of the criti- 
cal temperature with the applied flux are clearly visi- 
ble, which we associate with Little-Parks oscillations [s^l 
as already observed in many studies on superconducting 
wire networks and thin films with pinning arrays before 

We also plot the critical temperature vs magnetic field 
of a plain reference sample in Fig. 0] (a) and calculate 
from these data the coherence length ^(0) = 16 nm of 
the niobium by fitting it to the bulk expression Bc2 = 
0o/[2^C(r)2] with^(r) = ^(0)(l-T/T,)-i/2 m. As the 
coherence length is significantly smaller than the BCS co- 
herence length for niobium = 39 nm our films are in the 
"dirty limit" with a free mean path of I = 1.37^(0)^/^0 = 
9 nm. Assuming that the coherence length is not changed 
by the antidot patterning process and remembering, that 
the minimal width of superconducting material between 
two holes is W ^ Dg — Da = 400 nm, we calculate 
the reduced temperature, at which ^(IV) = W, to 
Tw/T^o = 0.9984. 

Above Tyy, which is marked with a dashed horizontal 
line in Fig.|4l the superconductivity in our system can be 
viewed as one-dimensional and the thin film expression 
= VT2(j}o/[2nWC{T)] should apply [HI, HI. This 
gives a parabolic dependence Tc/Tco oc for T/Tco > 
0.9984, what seems to be in reasonable agreement with 
our data. Tc/Tco of the perforated sample indeed appears 
like a Little-Parks modulated nonlinear background close 
to Tco with a development to a nearly linear behaviour 
for smaller temperatures T/Tco < 0.9984. 

By fitting the three single points at B = 0, B = Bi and 
B = —Bi to the above parabolic expression and using 
W = 400 nm, we calculate a coherence length ^(0) « 
18.2 nm in good agreement with the number extracted 
from the plain sample, although it is probably somewhat 
overestimated. The effective remaining superconductor 
width between the holes W is certainly larger than the 
used minimal value of 400 nm for two reasons. First, the 
holes have a circular shape and second, some of them are 
missing. So it might even be, that the coherence length 
in the perforated sample is somewhat smaller than in 
the plain, what would be supported by the impression 
of a slightly smaller slope in the more linear regime. A 
fit of the data for the perforated sample in this region 
however is difficult, as there are still modulations due to 
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FIG. 4: (Color online) Reduced critical temperature Tc/Tco vs 
normalized magnetic flux density B/Bi of a plain supercon- 
ducting thin film (blue line) and a sample with a quasihexag- 
onal lattice of antidots (red circles). The horizontal dashed 
line represents the temperature, at which ^{Tw /Tco) = W, 
the dashed parabola is a fit to the three data points of the 
antidot sample at B — 0, B = Bi and B = —Bi. 



the antidots superimposed. 

In summary we have investigated the Ic{B) — Tc{B) 
phase boundary of superconducting Nb films, which were 
patterned with quasihcxagonal arrays of submicron sized 
antidots. This fabrication method is especially interest- 
ing for microwave devices as it does not change the sub- 
strate and superconductor properties significantly. Our 
experiments revealed signatures of both, order and dis- 
order in the pinning lattice, which was patterned by us- 
ing a monolayer of self-assembling polystyrene colloids 
as microlenses for optical lithography. Pronounced sharp 
maxima of the critical current in the first and second 
matching field indicate a high ordering of the pinning 
sites, whereas the shoulder-like structures at noninteger 
higher field values might be related to pinning lattice 
blemishes and disorder. We also observe Little-Parks os- 
cillations of the critical temperature and approximately 
identify the wire network to thin film transition in our 
samples. The coherence length of the perforated film, 
extracted from the critical temperature vs magnetic field 
dependence, is in good agreement with that of a plain 
film, confirming that our patterning method has hardly 
influenced the properties of the niobium. 

We have performed experiments close to Tc, but when 
using superconductors with a higher magnetic penetra- 
tion depth, the results are also relevant at temperatures 
T = 4.2 K or even in the mK regime. This situation 
is given in very thin or dirty superconducting films and 
in different superconducting materials such as YBCO or 
NbN. In principle our patterning technique can be used 
with even smaller spheres, which would lead to commen- 
surability effects at much lower temperatures and higher 
magnetic fields. In these cases, the critical current and 
the pinning efficiency will modulate with the applied field 
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similar to the presented manner, what has to be consid- 
ered for the design of possible devices. 
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